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Abstract. Extended CCD monitoring of RR Gem revealed that it is a Blazhko type RRab star with the shortest 
Blazhko period (7?23) and smallest modulation amplitude (AM max < 0.1 mag) currently known. The short 
period of the modulation cycle enabled us to obtain complete phase coverage of the pulsation at each phase of 
the modulation. This is the first multicolour observation of a Blazhko star which is extended enough to define 
accurate mean magnitudes and colours of the variable at different Blazhko phases. Small, but real, changes in the 
intensity mean colours at different Blazhko phases have been detected. The Fourier analysis of the light curves 
shows that, in spite of the mmag and smaller order of the amplitudes, the triplet structure is noticeable up to 
about the 14th harmonic. The modulation is concentrated to a very narrow, 0.2 phase range of the pulsation, 
centred on the supposed onset of the H emission during rising light. These observational results raise further 
complications for theoretical explanation of the long known but poorly understood Blazhko phenomenon. 
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1. Introduction background of the phenomenon. Therefore we decided to 

begin an intensive photometric investigation of the star. 
Most RR Lyrae stars repeat their light curves with re- RR Gem (a = 07 ft 21 m 33!53, S = +30°52'59'/ 5, J2000) 

markable regularity. About 30% of the known galactic was discovered by L. P. Ceraski l)Ceraskilll90l| on pho- 

RR Lyrae stars, however, display cyclic modulation in tographic plates taken by Blazhko. The first visual ob- 

the shape and amplitude of the light curve over tens to ser vations of the star was made by iGrafJ lll905l) and he 

thousands of pulsational cycles (|Smith| [l995|) . The phe- was ab i e to determine the pulsation period (0^397). Since 

nomenon is often called the Blazhko effect since |Blazhko| tnat time a great num b er of visual observations have been 

CLMZl) was the first to recognize that the light maximum published, but there are onl y a few published photoelectric 

of an RR Lyrae star (namely RW Dra) showed phase mod- obser vations available (e.g.. lFitch et al.ll966tlLhi fc Janesl 

ulation on a long time-scale (around 42 days). Il989|) . 

In several well-studied RR Lyrae stars, the ampli- The absolute param eters of RR Gem were determined 

tude and the period of the Blazhko effect has been ob- bv iLiu fe Janesl l|l990j) from a Baade-Wesselink analysis. 

served to be variable on a time-scale of y ears l|Szeidllll988t It belon gs to the metal rich ([Fe/ Hl = -0.14 dex as com- 

ISmith et alJll999l lLaCluvze et al.1 12004). In this r espect piled in ljurcsik fc Kovacsl | l996) % l disk population group 

RR G eminorum provides a very interesting case. iDetrd of RR Lvrae stars l|Lavdenlfl995|l . 

lll970l) remarked that his unpublished photographic obser- Most of the papers dealing with Blazhko stars list 

vations from the 1930s demonstrated the Blazhko effect the possible theoretical explanations of the phenomena 

in RR Gem, but (in his words) it 'ceased to exist since and also the different observational arguments (see e.g. 

about 1940.' That the Blazhko effect might be temporary l LaCluvze et alJ l2004t lAlcock et all l2003l IChadid et all 

in nature may be crucial in understanding the physical |200J). Neither of the two groups of possible models 
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Fig. 1. Folded V light curve of RR Gem. The scatter of 
the light curve is caused by Blazhko modulation and not 
by observational inaccuracies. The amplitude of the mod- 
ulation is very small, it does not reach 0.1 mag even in 
maximum light, and the shape of the modulation points 
to a very regular behaviour. In this, and in the other fig- 
ures, phase 0.5 is set to the middle of the ascending branch 
defined as the phase where the V flux is equal to its time 
averaged value. 



a) the ma gnetic oblique rotator l|Shibahashil l2000t 
IConsenjlTflSaj) 

b) resonant cou pling of nonradial mode(s) to the do m- 
inant radial mode ijNowakowski k. Dziembowskill200l|) 
can fully explain the large variety of observational results 
of Blazhko behaviour. It is important to note, however, 
that both of the models connect the Blazhko period to 
the rotation of the star. 

In this first part of a series of two papers on RR Gem, 
we examine the photometric behaviour of the star exclu- 
sively using recent CCD observations. A reexamination of 
all the previous observations in order to reveal the long 
term behaviour of the modulation will be the subject of 
the second paper. 

2. The data 

The observations were obtained with the automated 60cm 
telescope of Konkoly Observatory, Svabhegy, Budapest 
equipped with a Wrigh t 750x1100 C CD (parameters and 
calibration are given in lBakosI <flQQSt> ^> using BV(RI) C fil- 
ters. The field of view was 17'x24'. Data reduction was 
performed using standard IRAF 1 packages. More than 
3000 frames were obtained in each passband on 56 nights 
between 14 January and 4 May in 2004 ( JD 2 453 019 - 
2 453130). The data were transformed to the standard 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Fig. 2. Amplitude spectrum of the V light curve (top 
panel). Middle panel shows the residual spectrum after 
prewhitcning with the pulsation frequency and its har- 
monics (fc/o, k = 1...15). The location of the modulation 
frequency (f m ) an d the symmetrically placed modulation 
frequency side lobes (fc/n ± f m ) are indicated by arrows. 
The spectral window is shown in the bottom panel. 



system by u s ing t he U,B,V,R,I magnitudes given by 
iLiu fc Janes! lll989h for SAO 060004 and SAO 060008. 
Relative aperture photometry to BD +31°1549 was calcu- 
lated, and SAO 060004 was used as check star to control 
the stability of the results. The magnitude differences be- 
tween our comparison and check stars remained constant 
with 0.0081, 0.0089, 0.0095, and 0.0093 mag r.m.s. scatter 
in the B, V, R, I colours, respectively. 

The magnitudes and colours we derived for 
BD +31°1549 (V = 10.277 mag, B — V = 0.366 mag, 
V - R = 0.233 mag, and V - I = 0.448 mag) are 
in excellent agreement with the V = 10.27 3 mag 
and B — V = 0.370 mag values given by ISturchl 



(1966). Photometric data are available elec- 
tronically at the CDS. In the electronic tables 
(http://cdsweb.u-strasbg.fr/cats/JA-l-A.htx ....) Column 
1 lists the HJD of the observation, and Column 2 gives 
the differential magnitude of RR Gem with respect to 
BD +31°1549 for the B, V, R, I colours. 

3. Photometric results 

The CCD observations revealed that, indeed, RR Gem 
shows Blazhko modulation, however, with very small, 
~ 0.1 mag modulation amplitude. The modulation pat- 
tern is very regular, with a definite node (fixed point) on 
the ascending branch. Fig. ^ shows the V band observa- 
tions superimposed on the mean light curve. The variation 
in the heights of maxima indicated a modulation period 
of about 7 days. 

The Fourier spectra of the V data and the data 
prewhitened with the pulsation frequency and its har- 
monics are shown in Fig. [21 Data analysis was performed 
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Table 1. Amplitudes and phases of the pulsation and modulation frequency components of RR Gem. 
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usin g the different applications of the MUFRAN pack- 
age l|KollatlJll990|) . The Fourier spectrum of the residual 
shows the triplet structure typical of amplitude modu- 
lation up to about the 14th harmonic component. The 
side lobe frequencies have similar amplitudes and are sym- 
metrically placed around the radial mode frequency and 



its harmonics. The amplitudes of the side lobes decrease 
much less steeply than the amplitudes of the harmonics of 
the main pulsation frequency (see Fig. . This fact may 
imply a different explanation of these frequency compo- 
nents than that of the radial mode where nonlinear effects 
invoke harmonic components with exponentially decreas- 
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Table 2. Side lobe amplitude ratios, phase differences, and their errors. Rk = j4fe/ +/ m/ Mfc/a— / m ! = 4>k} +f m — 
4>kf -f m - 
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Fig. 3. Amplitude ratios of the harmonic components 
of the dominant radial pulsation mode Ak/ a /Af (solid) 
compared to the amplitude ratio behaviour of the mod- 
ulation components A fc/o _/ m /A /o _ /m , A kfo+fm /Af 0+fm 
(open symbols) . While the amplitudes of the pulsation de- 
crease exponentially, the amplitude decrease of the modu- 
lation is linear. Error bars for the modulation components' 
amplitude ratios are indicated, the errors of the amplitude 
ratios of the dominant pulsation mode are smaller than the 
symbols used. 



This solution has been calculated by using locked val- 
ues of both the pulsation and modulation frequency com- 
ponents. The pulsation frequency adopted is the mean 
value obtained for the B, V, R, I data. The modulation 
frequency has been determined as the frequency which re- 
sults the smallest r.m.s. scatter of the residuals of the least 
squares solutions assuming symmetrical displacement of 
the side lobes up to the 9th harmonic component of the 
radial mode. The period of the modulation obtained this 
way is 7*?23. In Fig. 0] the V light curve is folded accord- 
ing to this modulation period; sinusoidal variation in the 
height of maximum light is evident. The ephemerides of 
maximum pulsation light and maximum pulsation ampli- 
tude determined from the measurements of the 2004 ob- 
serving season are the followings: 

HJD T max = 2 453 019.1590+0.3972884- E pulsation 

HJD T Blmax = 2 453 021.30 + 7.23 • E B i az hko. 

The error weighted mean value of the asymmetry pa- 
rameter (Q = ) averaged over the first 10 modula- 
tion component pairs in the four colours is —0.035 ±0.057, 
indicating very small asymmetry with slightly larger am- 
plitudes at the shorter frequency sides. For comparison the 
distribut ion of this paramet er of the MACHO data peaks 
at +0.3 l|Alcock et al.ll2003l) . 



ing amplitudes. The modulation components' amplitudes 
of AR Her (a Blazhko variable with large amplitude and 
phase modulation and a 32 d modulation period) show 
similar behaviour; from the fourth harmonic component 
its m odulation amplitu des dominate over the pulsation 
ones llSmith et al.lll999j) . 

Tabic n summarizes the amplitudes and phases of the 
pulsation and modulation frequencies in the B,V,R,I 
passbands calculated from the least squares solution using 
sine terms and HJD To = 2 453 19.0 initial epo ch value 
according to the formula given in iKovacsl l|l995l Eq. 2). 
In the top of Table ^ the mean magnitudes (An ) and the 
standard deviations of the fits are also given. 
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Fig. 4. The V light curve of RR Gem folded with the 7 d 23 
Blazhko period. The heights of maximum light clearly 
show a 0.1 mag amplitude modulation. 
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Fig. 5. V data folded with the Blazhko period (left panels) and their residual spectra after prewhitening with the 
pulsation period and its harmonics (right panels) for two equal length parts of the dataset. In the first half of the 
observations (between JD 2 453 019 and 2 453 174) the data coverage is nearly the same as for the whole data set and 
the residual spectrum is very similar to that shown in Fig. El for the complete data set (top panels). In the second 
part of the observations (between JD 2453077 and 2 453130) the data are more sparse, and the residual spectrum 
has changed significantly. The near equality of the sidelobe amplitudes has disappeared, instead, highly asymmetrical 
structure emerge, with very small if any power at the higher frequency sides (bottom panels). 



The currently viable theoretical explanations of the 
Blazhko phenomenon predict approximate equality of the 
amplitude ratios, and phase differences between the side 
lobe frequen cies of the different order components in 
each colour faovacslll995l IShibahashillgOOOt ISmith et all 
Il999t iNowakowski fc Dziembowskill200l|) . Although these 
predictions are based on simplified models it is worth 
checking the stability of these components in our data. 
Table El lists the amplitude ratios and phase differences 
and their errors of the side lobe frequency pairs in the 
B, V, R, I bandpasses for the first 10 modulation compo- 
nent pairs. Differences larger than 2a can be found in the 
amplitude ratios while the phase differences of the first 
5 modulation component pairs in the different colours 
do not agree within the 3er limit. The mean values and 
the standard deviations of the phase differences in the 
B,V,R,I colours of the first five modulation component 
pairs are 159.7±3.7, 144.7±3.5, 149.2±7.0, 148.7±5.5, and 
170.5 ± 2.3. The phase differences of the fifth components 
are statistically significantly larger in each colour than the 
phase differences of the second pairs. The large scatter of 
the third and fourth components indicates colour to colour 
changes of the phase differences, too. Differences between 
these pa rameters (Rk. Ad> k) have also been detected in 
AR Her ((Smith et alJll9 99). 



3.1. Stability of the results 

One of the theoretically unexplained characteristics of 
Blazhko stars is the high asymmetry of the amplitudes 
of the ± modulation frequency components. In contrast, 
only very small if any asymmetry can be seen in Fig. |2j 
We have checked what happens if data subsets with less 
complete phase coverage are used. The results are shown 
in Fig. [S] The 110 day dataset has been divided into two 
intervals of equal length but different number of observa- 
tions (~ 2000 and ~ 1000). In the first part of the obser- 
vations both the clear sky statistics were better, and the 
observational runs were longer. The residual spectrum of 
the first dataset is very similar to that of the complete 
dataset shown in the middle panel of Fig. El However, the 
residual spectrum of the second dataset is significantly dif- 
ferent, with enhanced power at the lower frequency side 
of the modulation, and with diminishing amplitudes at 
the higher frequency side. The data distribution of this 
second dataset also contains observations from each 0.1 
phase intervals of the Blazhko period, but for some phase 
intervals no complete light curve can be constructed. This 
is the typical data distribution of Blazhko star photom- 
etry, thus our result warns that classification schemes of 
the modulation based on the asymmetry of the side lobe 
amplitudes might be seriously biased by the phase cover- 
age of the data and may lead to an over-interpretation of 
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the results. As the Blazhko period of RR Gem is about 
an order of magnitude shorter than the typical period of 
the modulation, it should be expected that for stars with 
longer modulation periods even more extended observa- 
tions than we have on RR Gem are needed to correctly 
describ e the propert i es of the modulation. For compari- 
son, in lAlcock et all l)2003|) . 700 - 1000 measurements of 
the stars were analysed. 

The effect of data distribution biases on the side 
lobe amplitudes, however, has to be carefully treated. 
The fact that the different observations of RR Lyrae 
l|Szeidl k Kollatbl l2000t ISmith et alJl2003|) equally show 
the longer frequency modulation components to have 
higher amplitudes points to real differences in the am- 
plitudes of the modulation components. 

4. Variations during the Blazhko cycle 

The full coverage of the pulsation light curve at different 
phases of the Blazhko cycle makes it possible to study the 
changes of the light curves' shapes and Fourier parameters 
in detail, and to detect any possible variation in the mean 
magnitudes and colours. 



residual plot has been shown by iLee k. Schmidtl (j200lT) 
for DR And with the conclusion that 'large residuals are 
obtained from data points on the ascending branch'. To 
analyse such a behaviour with the Fourier decomposi- 
tion method might not be the most appropriate tool, as 
for a reliable model very high harmonic modulation fre- 
quencies are needed. This is reflected in the significantly 
different amplitude decrease behaviour of the modula- 
tion frequencies from the amplitude decrease of the har- 
monic components of the main pulsation (see Fig. 
Altho ugh no complete fully nonlinear mo del calculation 
exists. iNowakowski fc Dziembowskl l)2003|) showed that in 
evolved stars the nonradial modes are highly nonlinear. 
However, even in the case of high nonlincarity it is hard 
to explain why the high (5 — 10) order nonlinear coupling 
terms have amplitudes commensurable with the first order 
ones. 

It is also important to think about why the mod- 
ulation is concentrated in the minimum to maximum 
phase interval of the pulsation, where nonlinear effects are 
most important. In the pioneering spectroscopic study of 



4.1. Light curve changes 

We have divided the data into 10 subsets containing ob- 
servations from different phases of the Blazhko cycle. The 
minimum and maximum number of data in these subsets 
are 165 and 479, respectively. The behaviour of the Fourier 
parameters of the light curves in the different Blazhko 
phases is shown in Fig.H3 The total range of the variations 
in the A\, A\q amplitudes are 0.026, (between 0.390 and 
0.416 mag) and 0.009 mag (between 0.005 and 0.014 mag), 
respectively. The amplitude ratios (Rki = A^/Ai) fol- 
low the changes of the amplitudes but there is no signif- 
icant change in the epoch-independent phase differences 
(4>ki = 4>k — fc</>i) in accordance with the finding that no 
phase modulation in the light curve changes has been de- 
tected. 

We have checked how the light curves at different 
Blazhko phases deviate from theo mean light curve. Fig.0 
and [S] show the pulsation light curves at each one-tenth 
phase of the Blazhko cycle, and the residual curves after 
subtracting the mean light curve from the data subsets. 

Due to the good coverage of the data, a least squares fit 
with the pulsation period and its harmonics to the entire 
dataset (as shown in Fig. is a good representation of 
the mean light curve of the variable. The Fourier param- 
eters of this 'average' solution agree within the limits of 
the uncertainties with the corresponding amplitudes and 
phases of the complete solution listed in Table ^ 

The most striking property of the residual curves is 
that they show intense variations only in a narrow, ~ 0.20 
phase interval of the pulsation. The residual curve of all 
the data (Fig. ^ shows that the distortion is symmet- 
rical to phase 0.5 which is set to the phase of the mid- 
dle of the rising branch as described in Sect. 3. A similar 
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Fig. 6. Fourier parameters of the V light curves in the 
different phases of the Blazhko cycle. The amplitude ratios 
vary in phase with the A\ amplitude, but no changes in 
the epoch independent phase differences can be detected. 




Fig. 7. Pulsation light curve constructed from data falling 
into the 0.1 phase intervals of the Blazhko cycle. In each 
plot the mean light curve is also shown. Slight, but sys- 
tematic differences of the observations in most of the plots 
are evident. 



Fig. 8. The same data as shown in Fig. Rafter subtracting 
the mean light curve. Intense variations are concentrated 
in a 0.20 phase interval centred on phase 0.5 which corre- 
sponds in our representation to the middle of the ascend- 
ing branch, when H emission is supposed to occur. 



RR Lvrae. lPreston et al.l l|l965f) concluded that a displace- 
ment of the shock-forming layer in the atmosphere takes 
place during the Blazhko cycle. In IPreston fc Paczihskvl 
(1964) it was also shown that: 'the phases at which H 
emission occurs are closely related to phases of photomet- 
ric parameters'. Thus the phase of the U—B excess and the 
phase on the ascending branch when the visual luminos- 
ity of the star equals its time average value (at phase 0.5 
in the figures) are closely related, with only some minute 
differences, to the onset of the H emission. The U — B 
observations of iLiu fc Janea ((1989) indicate strong excess 
around the phase of the mid ascending branch in RR Gem. 
In Figs. |H1 and OH the symmetrical modulation is centred 
exactly on this phase of the pulsation, indicating a con- 
nection between the origin of the modulation and that of 
the H emission. 

As the standard deviations of the data from the fits 
using the 46 frequencies listed in Table ^ are larger in 
each colour than expected from the observations of the 
comparison and check stars, some further variations of 
the residual light curves are also suspected. The Fourier 
spectra of the residuals after prewhitening with both the 
pulsation and modulation frequencies listed in Tabled for 



the B, V, R, I data are shown in Fig. ^| Another modu- 
lation frequency component is present in each colour; the 
arrows indicate the positions of / = kfo + f' m (k — 1..6). 
This tertiary periodicity corresponds to modulation with a 
period of 6^69, which is very close to the dominant modu- 
lation period. A tertiary period shorter th an the dominant 
modu lation period has been also found bv lLaCluvze et all 
( 2004| for XZ Cygni. A quintuplet structure with equidis- 
tant frequency spacing would predict secondary modula- 
tion with a much shorter, 3^61 periodicity. 

Fig. El shows the residual V light curve after 
prewhitening with the frequencies of Table ^ folded with 
the pulsation period. Residuals are still concentrated in 
the same phase of the pulsation as in Fig. [§] This be- 
haviour of the residuals indicates that, perhaps, some ir- 
regular behaviour of the Blazhko modulation mimics the 
tertiary periodicity. This idea seems to be supported by 
separate analyses of the two halves of the data, as neither 
parts of the data show the same residual structure shown 
in Fig. [ill 
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Table 3. Magnitude and intensity mean brightnesses and colours in the different Blazhko phases. 









magnitude averages 






intensity averages 


Biphase Ai(V) 


B 


V 


R 1 B-VV-RV-7 


B 


V 


R 7 B-VV-RV-7 



0.0 - 0.1 0.3898 11.819 11.414 11.150 10.903 0.405 0.264 0.511 11.713 11.359 11.118 10.885 0.354 0.241 0.474 

0.1 - 0.2 0.3981 11.818 11.416 11.150 10.905 0.402 0.266 0.511 11.707 11.358 11.116 10.885 0.349 0.242 0.473 

0.2 - 0.3 0.4036 11.824 11.419 11.153 10.908 0.405 0.266 0.511 11.706 11.358 11.117 10.888 0.348 0.241 0.470 

0.3 - 0.4 0.4104 11.827 11.421 11.156 10.910 0.406 0.265 0.511 11.702 11.356 11.117 10.888 0.346 0.239 0.468 

0.4 - 0.5 0.4159 11.829 11.422 11.154 10.911 0.407 0.268 0.511 11.701 11.355 11.115 10.888 0.346 0.240 0.467 

0.5 - 0.6 0.4088 11.824 11.421 11.155 10.908 0.403 0.266 0.513 11.701 11.356 11.117 10.885 0.345 0.239 0.471 

0.6 - 0.7 0.4081 11.824 11.420 11.154 10.908 0.404 0.266 0.512 11.704 11.357 11.117 10.886 0.347 0.240 0.471 

0.7 - 0.8 0.4953 11.823 11.416 11.152 10.905 0.407 0.264 0.511 11.710 11.358 11.118 10.886 0.352 0.240 0.472 

0.8 - 0.9 0.3947 11.820 11.414 11.150 10.903 0.406 0.264 0.511 11.712 11.357 11.116 10.883 0.355 0.241 0.474 

0.9 - 1.0 0.3900 11.816 11.413 11.146 10.901 0.403 0.267 0.512 11.710 11.357 11.114 10.882 0.353 0.243 0.475 
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Fig. 9. The V residual light curve of RR Gem folded with 
the 0*?3972884 pulsation period after subtracting the mean 
pulsation light curve. 




frequency [c/d] 



Fig. 10. Amplitude spectra of the residual B, V, R, I light 
curves after prewhitening with the pulsation and modula- 
tion frequencies listed in Tabled Arrows indicate residual 
frequencies at fc/ + f' m (k = 1..6). 



Fig. 11. The V residual light curve of RR Gem folded 
with the 0^3972884 pulsation period after prewhitening 
with the pulsation and modulation frequencies listed in 
Table [T] Residual light variation is still concentrated in 
the same (0.4 — 0.6) phase interval of the pulsation as in 
Fig. El 

4.2. Magnitude and colour changes 

Multicolour photometry of Blazhko stars which were ex- 
tended enough to determine whether there is any change 
in the mean colour of the variables during the Blazhko 
period has never been obtained. 

According the mean characteristics of the Blazhko 
variables in the MACHO database, the modulation fre- 
quency occ urs in the spectra w ith an overall amplitude of 
0.006 mag l|Alcock et al.ll200"^) . i.e., this is the typical am- 
plitude of the the modulation in the magnitude average V 
light level of Blazhko variables. 

In Table|3]the magnitude and intensity mean B, V, R, I 
magnitudes of RR Gem in the different Blazhko phases are 
listed. Mean colours are calculated as the arithmetic dif- 
ferences between the mean magnitudes. Fig. 1121 shows the 
mean V magnitude and B — V and V — I colour variations 
during the Blazhko cycle calculated as magnitude and in- 
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tensity averages. Although the magnitude mean bright- 
ness level of the star seems to vary with a 0.005 mag semi- 
amplitude (in accordance with the MACHO results), the 
intensity mean brightness shows much smaller if any varia- 
tion. Concerning the colours, the situation is the opposite; 
the colours calculated from the intensity mean magnitudes 
definitely show correlated changes with the Blazhko phase, 
namely, the star is slightly cooler in its small amplitude 
phase than during the highest amplitude period. 

Whether intensity or magnitude mean colors better 
correspond t o the static value s of pulsating stars is still 
controversial ijBono et al . 1995). As the luminosities of the 
stars have real physical meaning, we think that the inten- 
sity means reflect the real parameters of the star more cor- 
rectly. According to the intensity mean values the mean 
brightness level of RR Gem remains the same during the 
Blazhko cycle, but a slight colour (temperature) variation 
takes place in the sense that the star is ~ 30 — 40K cooler 
when the amplitude of the pulsation is the smallest. 



5. Discussion 

The extended CCD observations of RR Gem exposed 
some previously unknown or unexamined properties of 
the Blazhko phenomenon. The new results concerning the 
properties of the modulation of RR Gem are summarized 
in the next items. 
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Fig. 12. Intensity and magnitude mean V magnitude and 
B — V and V — I colours of RR Gem in the differ- 
ent Blazhko phases. The intensity mean colours indicate 
colour variations during the Blazhko cycle, namely the 
star is slightly cooler in the smaller amplitude phase than 
during the highest amplitude period. 



First detection of a very low amplitude Blazhko modu- 
lation with symmetrically placed triplets in the Fourier 
spectrum. 



— The modulation is concentrated in a 0.2 phase interval 
of the pulsation, centred on the mid of the ascending 
branch. 

— In contrast to the exponential decrease in the ampli- 
tudes of the harmonic components of the radial pul- 
sation, the decrease in the amplitudes of modulating 
components is close to linear. 

— In contrast to theoretical expectations the amplitude 
ratios and phase differences of the ± modulation com- 
ponents are not the same for the different order mod- 
ulation component pairs and in the different colours. 

— According to the intensity mean colours and magni- 
tudes in the different phases of the Blazhko cycle, there 
is very small if any mean V brightness variation along 
the Blazhko period, while the mean colours indicate 
small but definite changes in the sense that the star is 
30-40 K cooler in its small amplitude phase. 

— The light curve cannot be completely fitted with sym- 
metrically placed Fourier frequency components alone. 
Residual scatter of the light curve is still concentrated 
in the ascending branch, indicating some irregular be- 
haviour of the modulation. 

Any model which relates the Blazhko periodicity to the 
rotation of the star predicts a similar percentage of small 
amplitude modulation as large amplitude ones. Because 
of the incidental angle of view of the rotational axes, 
the observed amplitudes of the modulation has to vary 
uniformly from zero up to some tenths of a magnitude. 
On the contrary, no small amplitude Blazhko modula- 
tion has been known previously. This fact has to be care- 
fully examined on the grounds of observational biases. 
The most extended observations of field RR Lyrae stars 
were the photographic observations from the middle of the 
20 th century. These observations were not, however, accu- 
rate enough to detect modulation amplitudes smaller than 
~ 0.2 mag. Photoelectric observations of field RR Lyrae 
stars are very sparse and usually time limited, so again, 
they were not suitable to detect small amplitude, long pe- 
riod modulations. The cluster and extragalactic variable 
surveys are usually seriously affected by crowding prob- 
lems and are also often very time limited, thus small am- 
plitude modulations might be missed or attributed to ob- 
servational inaccuracy. 

The detection of hundredths of magnitude amplitude 
modulation of RR Gem indicates that there might be a 
serious selection effect against small amplitude Blazhko 
stars. A systematic search using today's accurate pho- 
tometric detectors is needed to reveal the true unbiased 
statistics of the modulation amplitudes. 

Among the 731 fundamental mode Blazhko variables 
of the MACHO database there was only one star with a 
modulation per iod shorter than 1 days and only 5 with 
Pbi < 15 days (jAlcock et al.ll2003|) . In contrast, RR Gem 
is already the second very short modulation period vari- 
able in the much smaller sample of well observed galac- 
tic field Blazhko variables (less than 50 stars). The other 
short period Blazhko star (AH Cam) was studied in de- 
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tail bv lSmith et"afl l|l994() . Both AH Cam and RR Gem 
are relatively metal rich variables with pulsation periods 
shorter than 0^4. These are their common features, but 
the modulation pattern of the two stars are very dif- 
ferent. In AH Cam the amplitude of the modulation is 
0.3 — 0.4 mag, and phase modulation with about 0.1 pul- 
sation phases occurs as well. A simplified explanation of 
the larger percentage of short modulation period Blazhko 
stars in our Galaxy than in the LMC would relate the 
short modulation periods to the larger metal content of 
these st ars. However, the 530 day modulation period of 
RS Boo l|Kanvd|l986|) : lNagvllll998l)') wi th P puis = 0.377 d 
and [Fe/H]=— 0.35 dex IjLavderJ Il995[) indicates that a 
more compl ex description i s need ed. 

Though IChadid et ail l|2004h . favouring the nonra- 
dial mode explanation of the Blazhko phenomenon, ar- 
gue against the presence of a strong magnetic field in 
RR Lyrae, our detailed analysis of RR Gem revealed a 
challenge to these models as well. The two most impor- 
tant results which would be hard to explain by nonradial 
modes are the followings. 

a) The linear decrease of the amplitudes of the modu- 
lation components compared to the exponential decrease 
in the amplitudes of the radial mode (Fig. [3J indicates sig- 
nificant differences between these frequency components. 
As nonlinearity is proportional to the amplitude, it is ex- 
pected that if the harmonic terms of the radial mode in- 
voked by nonlinear coupling decrease exponentially, then 
the nonlinear coupling terms of the much smaller ampli- 
tude nonradial modes have to decrease even more steeply. 
The linear decrease of the amplitudes of the very small am- 
plitude modulation components indicates a different origin 
of its higher order components than nonlinear coupling. 
This anomalous behavior of the amplitudes of the modu- 
lation frequencies arises from the high concentration of the 
modulation in only 0.2 phase intervals of the pulsation. 

b) The modulation is centred on the phase when H 
emission is assumed to occur. Such a connection between 
a temporal physical phenomenon in the atmosphere and 
the occurrence of the modulation suggests that during the 
Blazhko period real changes in the atmosphere are present. 
This is also supported by the detected colour variations 
during the Blazhko cycle. 

Our detailed photometric observations of RR Gem 
have given new insight into the Blazhko phenomenon. 
However, it is an open question whether the Blazhko prop- 
erties of RR Gem are unique or give general information 
about the modulation as well. To find the answer, further 
similar studies of other Blazhko variables are planned. 
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